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CROSS REFERENCE TO RELATED APPLICATIONS 



This application claims priority under 35 U.S.C. §119 (e) to 
Provisional Patent Application Serial No. 60/424,141 filed 
10 November 6, 2002; the disclosure of which is incorporated by 
reference herein. 



This application is directed to the analysis of surfaces and 
in particular to the analysis of surfaces by calculating the 

20 conformal structure of the surface by providing a fundamental 
geometric tool for the analysis of surfaces by converting compact 
Riemann surface theory to computational algorithms . 

Geometric surface classification and identification are 
fundamental problems in the computer graphics and computer aided 

25 design fields. As scanning and imaging technology has developed, 
large numbers of colored meshes are becoming available in 
databases and on the world wide web (WWW) and the Internet. In 
addition, medical imaging technology, such as MRI and PET imaging 
systems are capable of producing three-dimensional (3-D) models of 

30 internal body structures. For example, recent developments in 
brain imaging have accelerated the collection and storage of such 
images in databases of brain maps. Similarly, in biometric 
security applications, face recognition involves the imaging, 
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storing, and matching of 3-D facial features to previously stored 
faces. Also, entertainment systems that use 3-D webpages are 
increasing in number, and computer animation techniques , such as 
morphing and texture mapping, also involve the creation and 
manipulation of 3-D surfaces. 

In all of these applications, the geometric data are 
represented as triangular meshes that' have a combinatorial 
structure instead of a differential structure. Accordingly, it is 
difficult to process these surfaces using differential geometry 
techniques. • Current analysis methods measure the Hausdorff 
distance between two surfaces; however, there is no general 
approach to find correspondence between the surfaces and in 
addition, combinatorial searching is inefficient. In addition, 
the current methods of surface analysis are heavily dependent upon 
the triangulation and resolution of the surface.. However, 
different triangulations and resolutions can result in widely 
varying results. Finally, geometric surface data are extremely 
large. One surface can have millions of vertices and faces such 
that the sheer number of calculations that are needed for current 
systems make it extremely difficult to develop effective and 
efficient algorithms. In addition, currently there is no 
effective general method to classify surfaces using topological 
invariants since the classification is typically too coarse, or 
using Euclidean geometric invariants since the classification is 
too rigid. 

Accordingly, it would be useful to provide a geometric 
analysis method that is an intrinsic system in that it depends 
only upon the geometry of the surface; that provides for a general 
way to classify surfaces effectively, find correspondence between 
two surfaces in the same class; and that provides for efficient 
and achievable computation that is both numerically stable and 
accurate . 
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BRIEF SUMMARY OF THE INVENTION 
A method for analyzing, classifying, and recognizing 
geometric surfaces is disclosed. Geometric surfaces are treated 
5 as Riemann manifolds and the conformal structure corresponding to 
the surfaces is calculated. The conformal structure of the 
surface contains the intrinsic geometric information about the 
surface, but in a much more compact format. In general however, 
surfaces are represented as a plurality of mesh data, with the 

10 number of mesh data points being quite large. Calculating the 
conformal structure of such a meshed surface can be a difficult 
undertaking due to the large number of mesh data points and the 
even larger number of calculations that are required. Conformally 
mapping the surface to a canonical parameter domain, such as a 

15 disk, sphere, or plane retains the geometric information of the 
surface, and renders the calculation of conformal structure much 
easier . 

In particular, in one embodiment, first and second surfaces 
are conformally mapped to a canonical parameter domain forming 

20 first and second mapped surfaces. The conformal parameterization 
for each mapped surface are computed and compared with one another 
to determine if the surfaces match. 

In another embodiment, a method for classifying a surface is 
disclosed in which the surface is classified according to the 

25 conformal parameterization. In particular, the period matrix R 
corresponding to the surface is determined and stored. 
Subsequently, a search for a particular surface can be conducted 
by examining the previously stored period matrix R and comparing 
this matrix to a second period matrix R' that corresponds to a 

30 desired surface. 

In another embodiment, a method for surface recognition is 
provided. In particular, a mesh representing a surface is 
provided and one or more feature points are sequentially removed. 
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For each feature point that is removed the corresponding period 
matrix R is calculated. By comparing the resulting sequence of 
period matrices to previously calculated sequences of period 
matrices, a surface may be recognized. 

Alternatively, all feature points can be removed at once and 
a point is selected within the surface. As this point is moved 
about the surface on a predetermined orbit, a sequence of period 
matrices are calculated and compared to a previously calculated 
sequence of period matrices. 

In another embodiment, a method of image compression is 
disclosed. A mesh representing a surface is provided and the 
conformal parameterization for the mesh is calculated. Using the 
conformal parameterization, the mean curvature can be calculated 
and with these two parameters, the original surface can be 

15 uniquely determined. 

In another embodiment, applications to medical imaging are 
disclosed. A medical image, such as of the brain or other organ 
is typically a genus-zero surface. Conformally mapping the genus- 
zero surface to a sphere enables the surface to be analyzed. 

In another embodiment, a method for animating a surface is 
disclosed. Given two similar shapes the feature points are 
removed from each surface and the doubling of each surface is 
computed. Each surface is decomposed to one or more patches and 
each patch is mapped to a plane. A conformal mapping from one 
25 plane to another is determined and after selecting control points, 
and a BSpline or other smooth curve function is used to generate a 
smooth transition between the two planes. 

In another embodiment, a method for generating textures to 
cover a given surface is provided. The surface is mapped using 
conformal parameterization to a canonical parameter surface, such 
as a plane surface, and the texture calculated for that parameter 
surface. To make the texture patches globally smooth, the 
Dirichlete method is used to diffuse the boundaries between 
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texture patches. In this way, the texture patches are * grown" and 
"stitched" together and then mapped to the parameter surface. 

Other forms, features, and aspects of the above-described 
methods and system are described in the detailed description that 
follows . 

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS 

The invention will be more fully understood from the 
following detailed description taken in conjunction with the 
accompanying drawings in which: 

Figs, la and lb depict a conformal mapping between a human 
face and a square; 

Figs, lc and Id depict a checker board texture mapped from 
the human face of Fig. la to the plane of Fig. lb; 

Figs. 2a-d depict various components of a holomorphic 1-form 
of a two hole torus; 

Fig. 3 is a spherical conformal embedding of a gargoyle 
model in a sphere; 

Fig. 4 depicts a human brain model conformally mapped to a 
sphere; 

Fig. 5 depicts a bunny model mapped to the unit sphere; 
Figs. 6a-b depict zero points of parameterization; 
Figs. 7a-d depict a global conformal atlas for genus two and 
three tori; 

Figs. 8a-d depict the topological equivalence but not 
conformal equivalence of two genus-one tori; 

Figs.9a-d depict genus-one surfaces with different conformal 
structures; 

Figs. 10 a-d depict the improvement in uniformity of the 
global conformal parameterization; 

Figs, lla-d depict various genus-two surfaces with different 
conformal structures; 
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Figs . 12a-b depict the use of regulariztion of the 
triangulation of a bunny surface; 

Figs. 12c-d depict a reconstruction of the bunny surface 
from a confonnal geometric image; 
5 Fig. 13a depicts a brain surface model; 

Fig. 13b depicts the brain surface model of Fig. 13a 
conformally mapped to a sphere; 

Fig. 13c depicts a spherical geometry image of the brain 
surface model of Fig. 13a; 
10 Fig. 13d depicts a brain surface model reconstructed after 

Fig. 13c has been compressed 256 times; 

Fig. 14 depicts a geometric morphing from a human female 
face to a human male face using conformal structures; 

Figs. 15a-b depict the global parameterization of a tea pot 
15 model at an original level of triangulation; 

Figs. 15c-d depict the global parameterization of a tea pot 
model at a simplified level of triangulation; and 

Figs. 16a-d depict the global parameterization results for 
four high genus surfaces. 

20 

DETAILED DESCRIPTION OF THE INVENTION 
In the embodiments that follow, two-dimensional (2-D) 
surfaces are treated as Riemannan surfacesand the conformal 

25 structure corresponding to the surfaces is calculated. All 
orientable surfaces are Riemann surfaces, and have an intrinsic 
conformal structure that is invariant under conformal 
transformations. . In general, the conformal structure is more 
refined than a topological structure and less rigid than a metric 

30 structure. For a genus-one surface, the space of all the 
conformal structure is two-dimensional. Thus, by using two 
parameters, all genus-one surfaces can be classified. In general, 
for a genus g surface, the space of all the possible conformal 
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structure is 6g-6 dimensional. Thus, all genus g surfaces can be 
classified using a g by g complex matrix. 

A methodology is provided to systematically compute the 
conformal equivalence between two surfaces is provided. In 
particular, for any two surfaces that have the same conformal 
structure, a method is provided to systematically compute the 
conformal one-to-one mapping between the two surfaces. For a 
genus-zero surface, the group of such mapping is 6-dimensional . 
For a genus-one surfaces, such a group is two-dimensional. For 
surfaces of more than genus-one, such a group of such mapping ■ 
includes only one dimension. Thus, advantageously, the methods 
described below provide an efficient method to find the best 
mapping and measure the Hausdorff distance between any two 
surfaces with the same conformal structure. 
15 For the methods described below in which the conformal 

structure of a surface is determined, the conformal structure is 
only a function of the geometry of the surface. It is unaffected 
by either triangulations and resolution and in addition, conformal 
mapping preserves the shape of the surface. 

It is well known that all surfaces are Riemann surfaces. 
Any Riemann surface has a conformal coordinate atlas, or a 
conformal structure. A conformal transformation maps a 

conformal structure to a conformal structure. Angles are 
preserved everywhere by a conformal transformation between two 
Riemann surfaces. As is known, a one -dimensional connected 
complex manifold is known as a Riemann surface. By Riemann 
uniformication theorem, all surfaces can be globally conformally 
embedded in a canonical space. The canonical space is typically 
a disk, a plane, or a sphere, the choice being determined by the 
intrinsic geometry of the surface. The conformally embedded 
surface includes a large portion of the original geometric 
information embedded onto the canonical spaces. Through 
conformal embedding, 3D surface matching problems can be 
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converted to 2D matching problems in these 3 canonical spaces. 
As discussed in more detail below, this method has the potential 
for non-rigid, deformed surface matching. 

The way of embedding the surface to the canonical space 
5 reflects the conformal structure of the surface. Specifically, 
all the global conformal embedding from a surface to the 
canonical space form a special group. If two surfaces can be 
conformally mapped to each other, they share the same group 
structure. In other words, such group structures are the 
10 complete conformal invariants. Hence, we can classify all 
surfaces using conformal invariants. For each topologically 
equivalent class, there are an infinite number of conformal 
equivalent classes. This is valuable for surface classification 
problems . 

15 Let Si and S 2 be two regular surfaces, parameterized by (x 1 , 

x 2 ) . Define a map $:S X -> S 2 represented in the local coordinates 

as <p (x 1 , x 2 ) = «39 1 (x l ,x 2 ) r (jPCxSx 2 )) . 

Let the first fundamental forms (Riemann metrics) of Si and 
S 2 be: 

20 dsf='£g u dx i dx J (1) 

u 

dsl^g (j dx i dx j . (2) 

U 

The pull back metric on Si induced by (p is 
mn a ax dx 

If there exists a function A(x l ,x 2 ) , such that 
25 ds\ =X{x\x 2 )fds 2 2 , (4) 

then we say that <f> is a conformal map between Si and S 2 . In 
particular, if the map from S z to the local coordinate plane (xi, 
x 2 ) is conformal, then (xi, x 2 ) is a conformal coordinate of Si, 
which is also referred to as an isothermal coordinate. Fig. la 
30 depicts a conformal mapping between a human face and a square on 
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the plane. Fig. lb depicts the conformal nature of the mapping by 
texture mapping a checkerboard to the surfaces. Inspection of 
Figs. la and lb illustrates that all right angles on the 
checkerboard are preserved on the texture of the surface. Fig. 16 
5 depicts the global parameterization results of four surfaces 
having a high genus, i.e., a surface with a genus >1. As can be 
seen, all angles on the checkerboard pattern are right angles, 
indicative of the conformal nature of the mapping. 

For a complex manifold, suppose U<z.C is an open set and let 
10 f be a complex function /:C/-»C. Then f is said to be 
holomorphlc, if for any z 0 g U there exists an e > 0 such that on 
the disk 

D(z 0 ,e) = {z e C\\z - * 0 1 < f } , ( 5 ) 

then f can be represented as a convergent power series 
15 /(z) = ]Ta,(z-z o y . (6) 

Let UczC and V <z Cbe open sets of C. A map f:U->V 

is biholomorphlc if f is one-to-one and holomorphic and f~ l :V->U 
is also holomorphic. 

Let S be a connected Hausdorff space with a family {[l/y,^)} 

20 that satisfies the following three conditions: 

1. Every Uj is an open subset of S, and S = {J JeJ Uj • 

2. Every Zj is a homeomorphism of Uj onto an open subset Dj 
in the complex plane. 

3. If Uj nU k & $ , the transition mapping 

25 Zv-ZvOzfiZjffJjnU^zfajnUt) (7) 

is a biholomorphic mapping, which is also a holomorphic 
homeomorphism. 

Thus, ^Pj> z jJij e j is ^ system of coordinate neighborhoods on S 

and . defines a one-dimensional complex structure on S. The 
30 coordinate neighborhood (U,z) of a Riemann surface is a pair of an 
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open set U in S and a homeomorphism z of U into the complex plane. 
U is referred to as a coordinate neighborhood of S and the 
homeomorphism z is referred to as a local coordinate or a local 
parameter. 

5 In general, a mapping f of S onto a Riemann surface R is 

said to be a holomorphic mapping, if wo/oz" 1 is holomorphic for 
all coordinate neighborhoods (U,z) of £ and (V r w) of R with 
f(u)czV . A biholomorphic mapping f:S-+R means that a 

holomorphic mapping f of S onto R has the holomorphic inverse 

10 mapping f l :R->S . 

Thus, two Riemann surfaces S and R are biholomorphic 
equivalent if there exists a biholomorphic mapping between them. 
If such a mapping exists, then S and R are regarded as the same 
Riemann surface and S and R have the same conformal structure. In 

15 general, complex structure, biholomorphic mappings and 
biholomorphic equivalence are also said to be conformal 
structures, conformal mappings and conformal equivalence, 
respectively . 

Let a surface S have a Riemann metric equal to 

20 ds 2 ^^^gijC&dx 3 , then the metric can be used to uniquely determine 

y 

a conformal structure {(Ui,z A )}, such that the local representation of 
ds 2 on a coordinate neighborhood (Ui,Zj.) is 

ds 2 = X{z i )dz i dz i , (8) 
where A(z { ) is a positive real function. 

25 To compute the conformal structure all the holomorphic 

differential forms on S must be found. Let S be a Riemann 
surface, then a holomorphic differential form CO on S is given by 
a family of { (U±, Zi,u>i) } that satisfies the following two 
conditions : 

30 1. Suppose {t/^z,)} is a conformal structure, then (a) has 

local representation OJi such that 
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where f± is a holomorphic function on Ui. 
2. If z,=^..(z y ) is a coordinate transition on £/,nt/ y *^ 



then 



/,W*,))#^/,<z,). 



A, — (10> 

Thus, the local representation satisfies the chain rule 
Mz^dz. =fj(zj)dzj . (11) 

The set of all holomorphic differentials on S is denoted as 
Q, 1 (S), where £l l (s) has a group structure that is isomorphic to the 
10 cohomology group of S. Thus, to compute Q}(s) the homology group 
of S must be computed. 

Let S be a two-dimensional Riemann manifold with metric 
g, N<=R 3 a compact two-dimensional manifold. For a C 1 map 
# = fa l ,f 2 ,f):S->.NcR* , let 

be the energy density in local coordinates x=(x 1 ,x 2 ) on S, 

S^{g a fi\iS afi ) = (gafi)' X • A C 1 variation of is a family (^) of C 1 map 

<f> e :S-±N smoothly depending on a parameter |fi]<So' an <* such that 

^ o =0. A variation fa) of is compactly supported if there 

20 exists a compact set Qc5 such that 0 e =0 on S/co on all |^|<^ 0 . 

A harmonic map on C 1 is a map 0:S-+Nc:R* that is stationary 
for Dirichlet's energy with respect to compactly supported 
variations and is given by 

J«M**j • d3) 

25 In local coordinates dA s = ^\g\dx l dx 2 , where |s|=det(g^). A map cp is 
harmonic if and only if 

A S t=An°0 r (14) 
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where X is a function globally defined on S and no$ is the normal 
at the image point on N. For a genus- zero surface, a harmonic 
mapping is a conformal mapping. If AT is R then cp is called a 
harmonic function. Note that all conformal maps are harmonic, but 
5 not all harmonic maps are conformal. 

A real differential 1-form t on S is harmonic if for any 
point on S, there exists an open set DaS , such that 
r\ D =df\ D , (15) 
where f is a harmonic function on S and d is the exterior- 
10 differential operator. 

All harmonic differentials form a special group H that is 
isomorphic to the cohomology group ^(SyR). According to Hodge 
theory, in each cohomology class, there is a unique harmonic 
differential form. 

15 A holomorphic 1-form o> can be decomposed into two real 

differential 1-forms x and y, such that a> = r , and in which 

both t and y are harmonic. By integrating a holomorphic 1-form 
on the surface, the surface can be conformally mapped to the 
complex plane. 

20 All holomorphic 1-forms form a group Cl x (S) that is the dual 

to the homology group Hi(S,Z). For a genus g surface S, there are 
2g generators of Hi(S, Z). Corresponding to each handle, there are 
two generators y i9 y {+g such that 

YrYi± z = 1,2, ...,g, (16) 

25 where " represents the algebraic intersection number of two closed 
curves. Then {^i^v-T^s-i^*} is called the canonical homology 
basis. If 6=^,^2,...^^,^} is a basis of H X (S,Z), the dual 
holomorphic 1-form basis is B* = {a i9 & 2 ...,& 2g _ l9 a> 2g } , satisfying 
Re faj =-r r i'j . (17) 
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Figs. 2a-2d depict the homology basis of a two-hole torus in 
Fig, 2a which consists of four closed curves. Fig. 2b depicts the 
harmonic 1-form co dual to e x in which the shaded curves are the 
integration lines of (O . Fig. 2c depicts the conjugate harmonic 
5 1-form co* that is orthogonal to harmonic 1-form depicted in Fig. 

2b. Fig. 2d depicts the holomorphic 1-form co + ^—lco* . 

The complete invariant for conformal equivalence is provided 
by a complex matrix. Suppose 2J = {p »72*} *- s a canonical 

homology basis and B* = {a> x ,6) 2 ...,G) 2g ^,co 2g } is a basis of Cl x (S), then 
10 the matrix P=(p±j) is called the period matrix of S, where 

Yi 

Examination of the period matrices of two surfaces , given by Pi 
and P2r respectively, can determine whether the two surfaces are 
conformally equivalent to one another, without the need to compute 
15 the conformal mappings between the two surfaces. 

In general,, surfaces are represented by triangular meshes. 
Every simplicial surface has a natural underlying complex 

structure. Let K be a simplicial complex, and a mapping / :|,Kj->2? 3 
embeds \K\ in R 3 , then M=(K, f) is called a triangular mesh, and K n 
20 where n=0,l,2 are the sets of n-simplicies . a* denotes the n- 
simplex, cr n = {v l9 u 29 ..JJ n } , where o t eK 0 . 

A chain space is the linear combination of simplicies and is 
given by 



C„(M) = 



2c y cr;|c # 6Z,cr;cJ:.J. (19) 



25 The elements in C nf n=0,l,2 are called an n-chain. Also, the 

summation of all faces y £ tk fk ^ s i n c 2 and M is also used to denote 
this 2-chain. 

A boundary operator d n : C n -> C B-1 among chain spaces is a 

linear operator. Let creK n9 a -[o Q9 o l9 ... 9 o n ^], then 
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d n <r n =2(-l)'[u 0i ... i u M> u Wf ..-u Il . l J . (20) 



Then for an n-chain in Cn, the boundary operator is defined as 

5 To denote the null space of d lr kerd 1 cz C t represents all the 

closed loops on M. Similarly, imgd 2 c:C x represents the image 
space of c?2 representing all the surface patch boundaries. Since 
di ' dz = 0, then imgdz <= kerdx. Hence, the homolgy group of M, 
Hn(M,Z) is given as 

10 H n (M 9 Z)= kera " . (23) 

Hi (M, Z) represents all the closed loops that are not the 
boundaries of any surface patch on M. The topology of M is 
determined by Hi (M, Z) . 

Let M be a closed mesh of genus g f and B = fa l9 y 29 ... 9 y 2g } be an 

15 arbitrary basis of its homology group. Then the intersection 
matrix C of B is given by 

c g = -/, '7j (24) 
where the " denotes the number of intersections, counting +1 when 
the direction of the cross product of the tangent vectors of ei 
20 and ej at the intersectin point is consistent with the normal 
direction and -1 otherwise. 

A co-chain space is the set of home omorph isms between chain 
spaces to R and are given by 

C n (M) = Hom(C n 9 R\n = 0,1,2 (25) 

25 where Hom(C n , R) represents the set of all homeomorphisms between 
Cn to R. The elements of are called n-cochains or n-forms. A 

coboundary operator is defined as 6 n :C n — >C n+1 . Let a^eCbe an n- 
form and c n+ ie C^+i is an n+l chain, then 

{fn&n X*« + l ) = «>n fo + ,S„ + , )r (26) 

30 and 81 ' 8o=0. 
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The cohomology group if (M,i?) is defined as 
H n (M 9 R)= k&XSn . (27) 

1-forms in kerS 1 are called closed 1-forms and 1-forms in imgS 0 are 
called exact 1-forms. Two closed 1-forms are called cohomologous 
5 if they differ by an exact 1-form. Cohomology group H 1 (M, R) is 
isomorphic to homology group Hi (M, Z) . 

Integration of an n-form along an n-chain is defined when 
Cn e Cn and co n eC n , as 



10 The boundary and coboundary operators are related by the Stokes 
formulae 



A bilinear operator star wedge product a*:C 1 xC 1 -^C 2 is defined 
similarly. Let feK2, the lengths of three edges as 1 0 , li, l2r and 
the area of f as A, then 



«*>*>c a >=& n {c n ). 



(28) 




(31) 



20 



where 



Q = (©(e 0 ),^(e 1 ),^(e 2 )) 

r=(r(^o),r(^i)>H^)) 



(32) 



(33) 



and the quadratic form G has the form 



. f -4/ 0 2 ll +11-11 ll+ll-lV 
TJ" A 2 1? +11-11 • (34) 

S U 2 +'o 2 -A 2 / 2 2 +A 2 -/o 2 -4/z 2 J 



25 



The harmonic energy CD of a closed 1-form is given by 
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E{g>)= *Zw e a>(ef , (35) 
where 

w e = — (cotor + cot/?) (36) 
2 

and if e is a boundary edge, ee9 2 M, then e attaches to one face fo 
5 and then vr e is given by 

w e =icota . (37) 

A closed 1-form is called a harmonic 1-form if it minimizes 
the harmonic energy, that is if the Laplacian operator defined as 

A*>(m)= E^kvpd^v]) (38) 

10 is equal to zero. Thus, a closed 1-form is harmonic if and only 
if its Laplacian is zero. Let M have a homology basis \ri,r 29 ... 9 r 2g } 

and a harmonic 1-form basis \co l9 a) 29 ... 9 a) 2g } , if 

< r, , a) j >= -Yi ' Yj > U J = l ? 2 v ..,2g (39) 

where -y^y. is the algebraic intersection number of yi and yj then 

15 the homology basis and harmonic 1-form basis are said to be dual 
to each other. 

Let M be a mesh and N is a smooth surface in R 3 . A piecewise 
linear map u:M ->N c: R 3 maps all the vertices of M to N as 
u(K 0 )czN. (40) 
20 The harmonic energy of u=(u 1 ,u 2 , u 3 ) is given as 

£(*) = £45"J (41) 

a 

where E(8\ia) is the harmonic energy defined for the 1-form (8ua) . 
If u minimizes the harmonic energy, E(u) , then u is a harmonic map 
and satisfies the following condition 
25 Au = (A8u\ASu 2 ,ASu 3 )=Ano U , (42) 
where n is the normal field on N. 
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Given a harmonic 1-form ©, there is a unique conjugate 
harmonic 1-form co*. A holomorophic 1-form is defined as 

gj + V-F*/. (43) 
All holomorphic 1-f orms form a group Q 1 (M) that is isomorphic to 
H 1 (M,R) . The basis of Q 1 (M) can be constructed directly from a 
basis of the harmonic 1-form group. Given a harmonic 1-form group 
having a basis of fa l9 a> 2 ,... 9 a> 2g } , then the basis of Q 1 (M) is given by 

+ V-T&> * , co z + 4^la)l 9 ... 9 a> 2g 

Given ^{ruVi^^g-x^g} is a basis of H X (M,Z), and 

B* = {a> l9 a> 2 ... ^g-i'&ig] is the dual basis of Q. 1 (M) , then a matrix 
C2gx2g= (Cij) and a matrix S2gx2g=(sij) are defined as 

Cy =<Yi>®} > (44) 

s 9 =<rt>&j>* (45) 

Then the period matrix R of M is defined as 

CR=SI. (46) 

Where R satisfies E 2 ^-!. The matrices (C,R) determine the 

conformal equivalent class of M. In particular, for any two 

surfaces Mi and M 2 with (l?i,Ci) and (R 2 ,C 2 ), respectively, then Mi 

and M 2 are conformal equivalent to one another if and only if 

there exists an integer matrix N such that 

N^R, N = R 2 ;N T C X N = C 2 . (47) 

The conformal structure of a mesh of genus g >0 is a family 
of {(UuZi.)} such that 

1. U± is simply connected and is formed by the faces of M, 

2. M c r\U x . 

3. zt is piecewise linear, and there exists a holomorphic 
1-form © such that bz± |ui « a>|ui. 

For a genus-zero mesh, there are no holomorphic 1-forms. In 
this case, the genus- zero surface can be conformally mapped to the 
surface of the unit sphere S 2 and the conformal structure of S 2 
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can be used to define the conformal structure of M. Thus, a 
discrete harmonic map u:M->S z defines the conformal structure of 
M. For any surface, by cutting M along c e Ci, a topological disk 
Dm can be formed and with it a special 1-chain. This cut along c 
5 is referred to as a locus or cut graph, and D M is a fundamental 
domain of M. The choice of c is not unique and accordingly, 
neither is the fundamental domain. 

A conformal map u:Dm->C can be found by using a holomorphic 
1-form a> + 4^l& eQ^Af) . A base point u 0 €Dm is selected and for 
10 any vertex osDm an arbitrary path yeCi(DM) is chosen, such that 
diY=o-o 0 , then 

u(u)=<fi>,^>+V = T<fl>\r > • (48) 

As discussed above, all genus-zero surfaces can be mapped to 
a sphere and therefore, all genus- zero surfaces are conformally 
15 equivalent. All conformal maps from S 2 to itself, form a six- 
dimensional Mobius transformation group. Using stereo-graphic 
projection to map the sphere to the complex plane, all Mobius 
transformations are of the form 

aZ * b *ad-bc = l 9 a 9 b 9 c 9 d gC . (49) 
cz + a 

20 However, to compute a conformal map to map a genus-zero surface to 
a sphere, extra constraints on the Mobius transformations are 
needed to * make the solution unique due to the form of the Mobius 
transformations . 

Another difficulty is that the image of the map is on S 2 and 

25 not in R 3 . Accordingly, when the map is updated, the image should 
be moved in the tangent space of S 2 and not in i* 3 . 

Having established the foregoing, several algorithms are 
provided below to compute the conformal structure described above. 
Applications in computer graphics, computer vision, and medical 

30 imaging fields are described. 
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In the algorithm that follows, Algorithm 1, the conformal 
maps between an arbitrary genus-zero surface and a sphere is 
calculated. First, the image mass center must be computed and is 
of the form 

5 mc(<f>) = \(fdA M . (50) 

For the discrete case, the following approximation may be used 

W*» ■ ^ j^Sw**^ (51) 

where A [UfV/ w3 is the area of face [u,v,w] . 

Algorithm 1 can now be used to compute conformal maps of 
10 genus-zero meshes to S 2 . 

Algorithm 1: Conformal Parameterization of Genus 0 Meshes 

Input: A closed genus-zero mesh M 

Output: A global conformal map <j>:M— >-S 2 

1. Compute the Gauss map, mapping M to S 2 
15 2. Compute the Laplacian at each vertex u of M, A<|>(u) . 

3. Project A(|)(u) to the tangent space of (j>(u)e S 2 . 

4 Update §{u) along the negative projected A<(>(u) . 

5. Compute the center of mass of A<|>(u), mc (<|>) , shift the center 

of mass to the center of S 2 , renormalize <j>(u) to be on S 2 . 
20 6. Repeat steps 2-5 for all vertices , until the projected 

Laplacian equals zero. 



Figs. 3, 4, and 5 depict spherical conformal mapping for 
three different genus-zero surfaces. In particular, Fig. * 3 
25 depicts a gargoyle model conformally mapped to S 2 , Fig. 4 depicts 
a brain model conformally mapped to S 2 , and Fig. 5 depicts a bunny 
model conformally mapped to S 2 . 

With regard to computing the conformal maps between any two 
topological disks, all such mappings form a three-dimensional 
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group that is a subgroup of the Mobius group discussed above and 
is represented by 

<fi(z)=^-^ 9 aa-bb=l,a,beC . (52) 
bz + a 

In order to compute the conformal maps between a topological disk 
and a unit disk, a technique referred to as doubling is used. 

Doubling converts surfaces with boundaries to closed 
symmetric surfaces. Given a surface M with a boundary 3M, a 
symmetric closed face M is constructed such that M covers M 
twice. That is, there exists an isometric projection n : M— >M 
that maps a face /eM isometrically to a face /eM . For each 
face f EiM there are two preimages in M . Algorithm 2 computes 
the doubling of a general mesh M. 

Algorithm 2 . Compute Doubling of an Open Mesh 

Input: A mesh M with boundaries. 

Output: The doubling of M, M. 

1. Make a copy of M, denoted as -M. 

2. Reverse the orientation of -M. 

3. For any boundary vertex u € 3M, there exists a unique 
corresponding boundary vertex -u e 3 -M, and for any edge on ee SM 
there exists a unique boundary edge -ee d -M. Find all the 
corresponding vertices and edges. 

4 . Glue M and -M such that the corresponding vertices and edges 
are identical. The resulting mesh is the doubling M . 

Using the doubling technique described in Algorithm 2, the 
conformal mapping of a topological disk to S 2 can be directly 
computed. Since the doubling surface is symmetric, M and -M will 
be mapped to a separate hemisphere and using stereographic 
projection n a hemisphere of the sphere can be mapped to the unit 
disk. In this manner, a conformal mapping is computed that maps 

-20- 



WO 2004/044689 




CTYUS2003/035395 



between the topological disk and the unit disk D 2 . By applying 
the Mdbius transformation in equation (52), all possible conformal 
mappings may be computed. 

5 

Algorithm 3. Compute a Global Conformal Map from a 
Topological Disk to D 2 . 

Input: A topological disk M. 
10 Output: A global conformal map <j> from M to the unit disk D 2 . 

1. Compute the doubling M of M. 

2. Compute a global conformal map (j>:M->S 2 , preserving symmetry . 

3. Rotate §(M ) such that ty(d M) is the equator, 

4. Use stereo-graphic projection % to map the upper hemisphere 
15 to the unit disk. 

5. Output 7t°<|). 

For surfaces with non-zero genus , the holomorphic 1-form 
group Q 1 (M) , which is determiend by the topology of the surface, 
20 is important in computing global conformal parameterization for 
these surfaces. To compute this group, the homology basis is 
computed first, the dual harmonic 1-form basis is computed next, 
and then the harmonic 1-form is converted into a base holomorphic 
1-form. 

25 

Algebraic algorithms for computing homology and harmonic 1- 
forms are introduced. Given a mesh M, the corresponding homology 
basis is computed using an algebraic topology method. Let of eJT w 
and cr^" 1 e K n ^ , then define 
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+ 1 +<rr'e dcr? 
-1 -a^edc? 
0 ±<x n k - l £d<T? 



(53) 



Then the n-dimensional boundary matrix is defined as 

dn=&,<*rh (54) 
The homology basis is then formed from the eigenvectors 
5 corresponding to zero eigenvalues of the following operators 

D = d*d t +d 2 dl . (55) 

Algorithm 4. Computing Homology Basis for Mesh M 

Input: Mesh M. 
10 Output: Homology basis {ft^v^ftg} - 

1. Compute the boundary matrices for 3 if 5 2 - 

2. Compute the Smith normal form of the matrix D = d\d x +d z d^ . 

3. ' Find the eigenvectors of D corresponding to zero 
eigenvalues, to form {ft*^—*^*} • 

15 

All harmonic 1-forms form the cohomology group that is the 
dual of the homology group Hx{M,Z) . A harmonic 1-form is both 
closed and harmonic. According to Hodge theory all the harmonic 
1-forms form a linear space that is the dual space of the homology 
20 group- Also, each cohomology classhas a unique harmonic 1-form. 

Algorithm 5. Computing a set of harmonic 1-form basis. 

Input: A homology basis {ft, y 2 5...,ft^} of M. 

Output: A harmonic 1-form basis {p) l9 co 29 ... 9 co 2g } . 
25 1. Set the values of c{ = -ft -y J9 i 9 j = l,2,...,2g . 

2. Solve the following linear system for o&i 

ScD t =0 
Atf>. = 0 

-22- 



WO 2004/044689 




CT/US2003/035395 



3. Output {aj 19 a> 29 ...,G> 2g }. 

As an alternative to the algebraic approaches used above, 
the homology, cohomology, and harmonic 1-forms may be calculated 
5 using combinatorial algorithms as follows. 

Algorithm 6. Computing a fundamental domain of mesh M. 

Input: A mesh M. 

Output: A fundamental domain Dm of M. 

1. Choose an arbitrary face fo e M, let D^fo, d Dm = d f 0 , put 
all the neighboring faces of fd that share an edge with f 0 into a 
queue Q. 

2. While Q is not empty 

a. Remove the first face f in Q, let 

d f = eo + e x + e^. 

b. Dm = Dm U f . 

c. Find the first e± e df, such that -e± e d Dm, 
replace -e± in d by {e±+i, e i+ 2}, keeping that 
order. 

d. Put all neighboring faces that share an edge with 
f and not in Dm or Q into Q. 

3. Remove all adjacent oriented edges in d D M that are opposite 
to each other, i.e.,. remove all pairs {e k , -e k } from d Dm. 

25 The resulting fundamental domain Dm includes all faces of M 

that are sorted according to their insertion order. The non- 
oriented edges and vertices of the final boundary of Dm form a 
graph G that is referred to as the cut graph. 

For the cut graph, Algorithm 7 computes the corresponding 

30 homology generators that are also the homology basis of M. 

Algorithm 7. Computing a homology basis of M. 
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Input: A mesh M. 

Output: Homology basis {pi»7a»—>72*}- 



1. Compute the fundamental domain Dm of M and determine the 
corresponding cut graph G. 
5 2. Compute a spanning tree T of G, let G/T={e if ez, e2 g } . 
3; Choose a root vertex reT, depth first traverse T. 
4. Let 3 ei = ti-si, there are paths from root r to t± and s±, 
denoted as [r, t±] and [r,si] then connect them to a loop yi = 
[r, tj.]-[r, s±] , 

10 5. Output fan^'—'fzg} as a bas is of HxiGrZ) and Hi (M, Z) . 



20 1. Compute a fundamental domain Dm, and the cut graph G of mesh 
M and compute a spanning tree T, G/T={ ei, 02, e2 g } - 

2. Let a>i.(ei)=l and o>i(e)=0 for any edge eeT. 

3. Suppose that Dm is ordered in the way that DM={fi, f 2f fn} , 
reverse the order of Dm to {f n , fn-i,...,fi} • 

25 4. While Dm is not empty: 

a. retrieve the first face f of D M , remove f from Dm, d f - 
So + &i + e^- 



b. divide {e k } into two sets, F = {eedf\-eedD M } , 
Tl = {eedf\-e£dD M } . 




To explicitly compute a basis for the cohomology group of M, 
H 1 (M, Z) , a set of closed 1-forms \p) l9 a> 29 ... 9 <D 2g \ is found such that 

<Ya co j >=<?/ • (56) 
15 Where Sf is the Kronecker delta and y x is a homology basis. 
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c. choose the value of (D t {e k \e k eU arbitrarily, such that 
5] e6n <2>.(e)=-^ 6r ^(^)/ if II is empty, then the right hand side is 

equal to zero. 

d. Update the boundary of Dm, let dD M -dD M +df . 

5 Once the cohomology basis {a\ 9 a) 29 ... 9 G> 2g } has been computed, the dual 
of the homology basis, >/2* 

}, can be found by the linear 

trans form {a> l9 a> 29 ... 9 Q} 2g } such that 

10 Algorithm 9. Diffuse a closed 1-form to a harmonic 1-form. 

Input: A mesh M, a closed 1-form o. 

Output: A harmonic 1-form, cohomologous to ©. 

1. Choose f e C° (M) , such that A(<a + #0 = O. 

2. Solve the above sparse linear system for f. 
15 3. Output co+8f. 

Where 

A(»+<rX")= E^(4.v])+/(v)-/(4"^.- (58) 

[H,v]eM 

20 Given a harmonic 1-form {a^,0 2 ,...,0 2 *} , the conjugate harmonic 

1-form ©* can be found by solving the linear system 

2g 

<o>. ag> J9 M >=<o>, /\a>\M > . (59) 

Once the fundamental domain has been computed, the conformal 
mapping may be computed directly by integrating a holomorphic 1- 
25 form co. First, select a root vertex v 0 eDM then use the depth 
first search method to traverse the Dm. Each vertex usDm has a 
unique path y from v 0 to u, then we define 0(u)=<G>,y>. 
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Algorithm 10. Global Conformal Parameterization of a Mesh M 

Input: A mesh M f a holomorphic 1-form co. 

Output: A map 0:Dm-»C, or a global conformal parameterization. 
1. Compute a fundamental domain Dm of M. 
5 2. Use depth first search method to traverse the vertices ugDm, 
record the path from root vertex vo to u, denoted as Yu- 

3. Compute the integration 0(u) =«o, y u >. 

4. Output ^(u) as the conformal coordinates of u. 

10 Algorithm 11. Conformal Structure of a Mesh M 

Input: A mesh M. 

Output: A conformal structure of M { (Ul, z±) } 

1. Compute a holomorphic 1-form basis {ey, + V-TtfJ/* } . 

2. Compute a partition {U±} r such that M c U±, 01 is simply 
15 connected. 

3. For each U± choose a holomorphic base integrate 

the holomorphic 1-form on U±, denote the mapping as Zi. If there 
are zero points, subdivide U± and repeat step 3. 
4 . Output { ( U±, z±) } . 

20 

The global conformal parameterization obtained by 
integrating a holomorphic 1-form on a fundamental domain can be 
used for canonical decomposition of meshes, converting meshes to a 
tensor product spline surface, surface matching and recognition, 

25 and other useful image processing applications. 

According to Poincare-Hopf theory, a holomorphic 1-form a> 
must have zero points if M is not homeomorphic to a torus. Zero 
points of 03 are the points where the conformal factor is zero. A 
genus-g surface has 2g-2 zero points. A conformal mapping wraps 

30 the neighborhood of each point twice and double covers the 
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neighborhood of the image of p on the complex plane. Locally the 
map, ^: C— >C is similar in the neighborhood to 



Figs. 6a and 6b depicts the zero points on the global conformal 
parameterizations for an open teapot model and for the complex 
plane respectively . 

We can treat a harmonic 1-form co as a mapping from the 
surface M to the unit circle S 1 . Then for a holomorphic 1-form, 
the harmonic 1-form of the real part is the circle valued mapping. 
The harmonic 1-form of the imaginary part is the gradient field. 
The integration curves through the zero points will subdivide the 
surface into regular patches. In particular, for a mesh M and a 

holomorphic 1-form <2> = r + V-Ir* , the integration curve along x or x* 
and through the zero points partitions the surface into 
topological disks or cylinders. 

Let M be a topological torus M that is conformally mapped to 
C. By integrating a holomorphic 1-form co on its universal 
covering space, a periodic conformal map results. Selecting a. base 
point Uo, the image set of the base point is 
{a <A i ,a>>+b<Z 11 ,a>>+z 0 \ a,beZ). (61) 
This mapping is periodic, or modular. The entire torus is mapped 
into one period, which is a parallelogram spanned by <yi, co>, <yz, 
cd>, which are referred to as the periods of M. If the genus-g of 
M is greater than one, different handles may have different 
periods. The entire surface is mapped to g overlapping modular 
parallelograms. The parallelograms may attach to and cross each 
other through the image of the zero points. 

Figs. 7a-d depict this phenomena. In Figs. 7a and 7b a two- 
hole torus is separated into two handles and each handle is 



<f>(z)=z 2 



(60) 
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conformally mapped to a modular space. Figs. 7c and 7d depict a 
genus-three torus and the conformal mapping into modular space. 

To generalize the methods described herein meshes with 
boundaries will now be considered. Given a mesh M with 
boundaries, the doubling M of M is computed. For each interior 
vertex u e M, there are two copies of u in M f which are denoted 
as uj and u 2 . Ui and Uz are dual to each other as 
u x =w 2 , u 2 =Wj . (62) 

For each boundary vertex u e dM, there is only one copy in M, so 
that u is dual to itself. 

To compute the harmonic 1-forms on M, it is known that all 
symmetric harmonic 1-forms of M are also harmonic 1-forms on M. 
Define the dual operator for each harmonic 1-form a* as 

m([u 9 v])=co([u 9 v]). (63) 
Any co can be decomposed into a symmetric part and an asymmetric 
part as 

co = + (a -W) , ( 64 ) 

2 2 

where ±(a) + co) is the symmetric part and ^{co-a?) is the asymmetric 
part. 

Algorithm 12. Computing a set of holomorphic 1-form basis for 
meshes with boundaries 

Input: Mesh M with boundaries. 

Output: Holomorphic 1-form basis for mesh M of the form 

1. Compute the doubling of M, M. 

2. Compute the harmonic 1-form basis of M {a) l9 a) 2 ,...,(D 2g } . 

3. Assign Ti = + remove redundant ones. 

4. Compute conjugate harmonic 1-forms of x± denoted as t±*. 

5. Output the holomorphic basis \r t + 4-^t[ 9 t 2 + V-lT*,...,r^ -h-v/^Tr^}. 
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Figs. 8a and 8c depict two genus-one surfaces, that although 
they are topologically equivalent, i.e., both genus-one surfaces, 
the two surfaces are not conformally equivalent. Each torus can 
5 be cut open and conformally mapped to a planar parallelogram as 
depicted in Figs. 8b and 8d respectively. The shape of the 
respective parallelogram indicates the conformal equivalent class. 
The conformal equivalent classes are determined by the acute angle 
of the parallelogram, a right angle in these two cases, and length 

10 ratio between the two adjacent edges to represent the conformal 
invariants, or shape factors of these two genus-one surfaces. As 
depicted in Figs. 8b and 8d, the two tori have different shape 
factors and are not conformally equivalent. 

Table 1 below contains the conformal invariants of the 

15 genus-one surfaces depicted in Figs. 9a-9d. It is clear that none 
of the surfaces depicted in Figs. 9a-9d are conformally 
equivalent . 



Table 1: Conformal Equivalents of Genus-one Surfaces 



Mesh 


Shape Factor 


Vertices 


Faces 


torus 


1.0-1.1421 


1089 


2048 


knot 


1.0-.2721 


5808 


11616 


knot 2 


1.0+0.1281 


2050 


3672 


rocker 


1.0-3.5091 


3750 


7500 


teapot 


1.0-.1121 


17024 


34048 



20 Algorithm 13. Verify whether Ml and M 2 are conformally equivalent 

Input: Two meshes Ml and M 2 . 

Output: Indicia of the conformal equivalence or not of Mi and M 2 
1. Compute period matrices (Ri, Ci) and (Rz, Cz) corresponding to 
M x and Mz respectively. 
25 2. Compute the Jordan normal form of Rx = FiAPi*" 1 and 
Ri - PiTzPz- 1 . 

3. If Ti & Tz return false. 
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4. Let N = PiP2~ l r return true if N is an invertible integer 
matrix and NCilf = Cz, otherwise return false. 

The conformal factor X(u,v) indicates the first fundamental 
5 form of the surface S. If X is a constant then the Gaussian 
curvature of the surface is zero. By selectively cutting on the 
surfaces, new boundaries are introduced, thus the conformal 
structure can be altered. In practice, it is helpful to improve 
the uniformity of the parameterization and in general these cuts 

10 are made on the regions of the surface having a high Gaussian 
curvature. Figs lOa-d depict the improvement in uniformity. In 
the spherical parameterization depicted in Fig. 10a, the ear part 
is highly under sampled. By introducing topology cuts at the ear 
tips, the parameterization becomes much more uniform. 

15 In general the stability of the computations is highly 

dependent on the quality of the triangulation. If all angles of 
the triangulation are acute angles, the computing algorithms are 
guaranteed to be stable and convergent. Fig. 15 depicts the 
global parameterization of a tea pot model at two different levels 

20 of surface model complexity. As can be seen in Figs. 15a-b, for 
the more complex original tea pot the global parameterization 
results in all angles being acute angles and in particular right 
angles. Figs. 15c-d depict the global conformal parameterization 
of the simplified tea pot model in which all the angles are acute 

25 angles and in particular right angles. In both cases, regardless 
of the complexity of the model, the computing algorithms are 
convergent and stable. The following algorithm approximates a 
triangulation with all acute angles. 

30 Algorithm 14. Triangulation of a surface with all acute angles 

Input: A mesh M 

Output: Remesh M with all acute angles 

1. Subdivide the mesh using loop subdivision method. 
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2. Simplify the mesh using Edge Collapse using minimum edge 
length criteria. 

3. Repeat steps 1 and 2 until all angles on M are acute. 

4 . Output the remeshed M. 

Surface Matching Based on Conformal Parameter and Mean Curve 
Matching 

If one surface can be deformed into another one without too 
much stretching, such as human expression or skin deformation, 
then the deformation can be accurately approximated by global 
conformal mapping. Since conformal parameterization depends on 
the first fundamental form of the surfaces, and in particular the 
conformal structure depends on the Riemann metric continuously, as 
long as the Riemannian metric tensor does not change too much, the 
conformal structures are similar. Thus, mapping two surfaces to a 
canonical parameter domain and matching the surfaces in the 
parameter domain allows 3-D matching problems to be solved more 
efficiently. 

By storing the conformal factor Mu,v) and normal n(u,v) on 
the parameter domain, the original surfaces can be reconstructed 
uniquely up to rotation and translation in R 3 . X(u,v) defines the 
first fundamental form and n(u,v) defines the third fundamental 
form and hence the second fundamental form, i.e., the embedding in 
R 3 can be computed. Thus, the surface can be constructed uniquely 
up to a Euclidean transformation. 

A more efficient method is to use the mean curvature on the 
conformal parameter domain. For any closed surface without 
boundaries, the surface is uniquely determined by the conformal 
factor Mu,v) and mean curvature H. For any open surface with a 
boundary, the surface is uniquely determined by the conformal 
factor Mu,v), the mean curvature H, and the second fundamental 
form on the boundary. 
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To match surfaces based on Gaussian curvature and mean 
curvature, the surfaces to be matched are embedded in a canonical 
parameter domain. For example, a human face can be mapped to a 
unit disk. The Gaussian curvature and mean curvature are computed 
5 using conformal parameterization. The level sets of Gaussian 
curvature and mean curvature are families of planar curves on the 
parameter domain. These level sets of curves are then used to 
match the surfaces. 

To match surfaces that contain special features, the feature 
10 points are first removed and the doublings of the surfaces are 
computed. Next, the homotopy type of the map are constrained to 
guarantee that the features in the first surface are matched to 
corresponding features in the second surface. The conformal 
structures are then computed to perform the matching as described 
15 above. For example, to match human faces, the features such as 
the eyes, tip of the nose, and the mouth are removed prior to 
computing the conformal structure. 

Surface Classif icaiton 

20 To classify surfaces to allow efficient databasing and 

searching, the conformal structure in the form of the period 
matrices for each surface are computed and stored. Figs. 11 a-d 
depict various genus-two surfaces. As can be seen below, none of 
the surfaces depicted in Figs, lla-d are conformally equivalent as 
25 the period matrices R are not equivalent. 

The two-hole torus of Fig. 11a includes 8 61 vertices and 
1536 faces and has a period matrix R that is 
f-1.475e-3 4.840e-4 4.501e-l 2.132e-2 ^ 
4.858e-4 -1.439e-3 2.132e-2 4.501e-l 
-2.260e + 0 1.090e-l 1.467e-3 -4.858e-4 
1.090e-l -2.260 + 0 -4.840e-4 1.439e-3 J 
The vase model depicted in Fig. lib has 1582 vertices and 2956 
30 faces and a period matrix R that is 
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fl.053e-3 -8.838e-6 4.479e-l 2.127e-2^ 
-1.080e-4 -1.031e-3 2.172e-2 4.042e-l 
2.309e + 0 1.241e-l 1.053e-3 -1.080e-4 
v -1.241e-l -2.564e + 0 8.851e-6 1.031e-3 J 
The flower model depicted in Fig. 11c has 5112 vertices and 10000 
faces and a period matrix R that is 
t 6.634e-3 -1.950e-3 2.861e-l -6.0766-2^ 
-1.909e-3 7.091e-3 -6.076e-2 2.497e-l 
-3.768e + 0 -9.111e-l -6.634e-3 1.909e-3 
^-9.111e-l -4.303e + 0 1.950e-3 -7.091e-3j 
The knotty bottle depicted in Fig. lid has 15000 vertices and 
30000 faces and a period matrix R that is 
f-1.911e-2 2.757e-3 5.617e-2 -1.001e-3^ 
1.213e-3 - 9.294e-2 -1.003e-3 5.699e-2 
-1.792e+l -4.829e-l 1.912e-2 - 6.224e-4 
^-4.817e-l -1.819e + l -3.355e-3 9.295e-2 ) 



(67) 



(68) 



Surface Recognition 

It is desired that surfaces can be recognized without having 
to be matched to one another directly. Modifying the conf ormal 
structure of the surface in a canonical way and computing the 
period matrices for each modification provides a sequence of 
period matrices that indicate the intrinsic geometric properties 
of the surface. 

For example, to recognize a human face, the feature points, 
such as the center of the left eye, the center of the right eye, 
the nose tip, and center of the mouth are removed. For each 
modification to the facial surface, the doubling of the surface 
and the period matrices are computed. By comparing sequences of 
period matrices, we can recognize a geometric surface, such as a 
face. 

Alternatively, all the major feature points are removed and 
another point is selected and moved within the surface and period 
matrices of the doubling are computed for each movement of the 
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selected point. For example to recognize a human face the points 
at the center of the eyes, the tip of the nose, and the center of 
the mouth are removed and another point on the face is moved along 
a prescribed orbit. At each step, the point at the current 
5 location is removed and the period matrix computed. A sequence of 
period matrices will be computed, one for each point along the 
prescribed orbit. It is these period matrices that are used to 
recognize the surface. 

10 Harmonic Spectrum Analysis 

Alternatively, the Laplacian operator, described above, has 
infinite eigen values and eigen functions. The spectrum of all 
the eigen values reflects much of the intrinsic geometry of the 
surface. In addition, the eigen functions can be used to 

15 reconstruct the surface. Rhe surface can be recognized using only 
the spectrum of the surface as the signature of the surface. For 
example, in the medical field, by analyzing the spectrum of the 
shape of internal organs some illnesses may be detected. 

The desired eigen values and eigen functions can be computed 

20 for a surface represented by a triangular mesh by finding the 
eigen values and eigen vectors of the Laplacian matrix. 

Compression of Surface Data using Hamonic Eigen Functions 

A genus-zero surface is conformally mapped to the unit 
25 sphere and the position vector of the surface is represented as a 
vector valued function defined on the sphere. The eigen-f unctions 
of the Laplacian operator on the sphere are the spherical 
harmonics that form a basis for the functional space of the 
sphere. The position vector is then decomposed with respect to 
30 the functional basis and the spectrum is obtained. By filtering 
out the high frequency components, the surface data is compressed. 
Through the use of a MSbius transformation, described above , a 
region can be "zoomed into" for further examination. For general 
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surfaces, conformally mapping the surface to a canonical shape in 
its conforraal equivalent class and decomposing the surface 
position vector using the eigen-functions of the Laplacian 
operator provides the desired functional basis from which the high 
5 frequency components can be removed prior to storage. 

Alternatively, the conformal factor and mean curvature 
defined on the conformal coordinates can be used to determine the 
surface uniquely to a Euclidean transformation. In this method 
the two functions defined on the plane, i.e., the conformal factor 
10 and mean curvature, are used to represent the surface. Thus, a 
savings of one-third the storage is realized. Further compression 
may be obtained by using the eigen-f unction technique described 
above or other known compression techniques. 

15 

Rexneshing and Hardware Design 

By using conformal structure, we can remesh the surface 
after it has been conformally mapped into the parameter domain. 
In this way the irregular connectivity can be changed to regular 

20 triangulation. In theory, the reconstructed normals are accurate. 
This will simplify the representation of geometric data, and 
simplify the graphic hardware architecture . Currently, for 
general graphics hardware, there are memory buffers for storing 
connectivity information. The communication of between CPU and 

25 the graphics card that is necessary to specify this connectivity 
information is extremely time consuming. If the connectivity of 
the data stored within the memory on the graphics card is regular, 
and the graphics card can predict it by itself, then no extra 
memory will be needed for connectivity information. Thus, 

30 reducing the necessary level of communication between the 
processor and the graphics card. With respect to the architecture 
of graphics cards, currently the pipeline for processing the 
geometry and the pipeline for processing surface texture are 
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separate. Using regular connectivity, geometry can be represented 
as texture also, and these two separate pipelines can be combined. 
In this way, the complexity of the graphics card architecture can 
be reduced. 

Also by remeshing, geometry images can be constructed and 
the image format can be used to represent the surface geometry. 
In this way, many image processing techniques that operate on the 
geometry, such as compression, multi-resolution, and filtering, 

among others, can be used. 

Fig. 12a depicts a bunny model having irregular connectivity 
of the original mesh. After remeshing using the conformal 
structure, as depicted in Fig. 12b the connectivity is very 
regular and the reconstructed normals are very accurate. The 
conformal geometry image is shown in Fig. 12c, and the 
reconstructed shape is depicted in Fig. 12d. 

Parametric Surface and Mesh Conversion 

In the CAGD field, parametric surfaces such as BSpline 
surfaces and Bezier surfaces are used frequently. In the 
manufacturing industry, a controller that makes use of these kinds 
of parametric surfaces often guides the processing machines. 
However, geometric data are often represented as triangle meshes. 
Current geometric data acquisition devices output geometric data 
as dense point clouds. It is easier to convert these scanned 
point clouds into meshes and therefore it is very important to 
convert parametric surfaces to and from meshes. Currently, there 
are no automatic methods to convert meshes to Spline surfaces. 

By using the conformal geometry techniques described herein, 
this problem can be solved. As discussed above, a global 
conformal parameterization of the surface is computed and the 
surface is decomposed to canonical patches using the integrati6n 
lines along the gradient field through zero points. Each 
canonical patch is mapped to a rectangle on the plane, and a 
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tensor product spline surface is constructed on it. The resulting 
parameterization can be made globally smooth by matching the 
control points on the boundary. Hence, it is convenient to 
convert the mesh to parametric surfaces with an arbitrary 
5 desirable continuity. In addition, this construction preserves 
accurate normal information. 

Numerical Computation on Surfaces 

Conformal structure is a good parameterization for computing 
10 covariant differentiations on surfaces. Covariant differentiation 
is intrinsic to the surface geometry, so the embedding in the 
Euclidean surface is irrelevant. Conformal structural analysis has 
potential to compute natural physical processes on deformable 
surfaces . 

15 By using conformal coordinates, the differential operators 

have a very simple format. For an example the Laplacian operator 
is 

^=£"kh|-- (69) 
a £ dx a dx a 

This technique allows for the easier solution of surface 
20 partial differential equations such as Navier-Stokes equations and 
Maxwells equations. Using the conformal structure described 
above, the Gaussian curvature of a surface is easily determined. 

Medical Imaging 

25 The conformal structures described above can also be applied 

in the medical imaging field, such as in brain mapping, brain 
registration, heart surface matching, and vessel surface analysis. 
For example, by mapping the brain surface to the unit sphere, it 
is convenient to compare two brains and match the features. By 

30 analyzing the geometric structures on the brain, it is easier to 
find changes to a brain over time and to find potential illnesses. 
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The conformal map from a brain surface to a sphere is 
independent of triangulation and resolution. The conformal 
mapping provides a nice canonical space for us to compare and 
register two brain surfaces. Since the brain surface is very 
5 complicated, it is very hard for other methods to trace the 
evolution of the vertex's flow. The methods described herein 
handle the complicated surface structures while maintaining 
accurate angle information. Since the brain is typically a genus- 
zero surface, Algorithm 1, described above, may be used to map the 
10 brain surface to the unity sphere. Figure 14 shows examples for a 
brain mapping. 

Animation 

Conformal geometry can also be applied to computer graphics 

15 animation. Using current data acquisition technology, 3D shapes 
of an actor can be scanned with different gestures and 
expressions. Using the conformal analysis techniques described 
above, these key gestures and expressions can be mapped to one 
another. By using spline interpolation techniques, smooth 

20 transitions between the gestures and expressions can be generated 
between them. Thus, arbitrary shapes can be animated, including 
soft shapes and deformable models, which are extremely difficult 
to animate using current methods. 

Suppose we are given two similar shapes. First the feature 

25 points are located, and then removed. The doubling of the 
surfaces is computed and the homotopy type of the mapping is 
determined. A holomorphic 1-form on each surface is selected, 
such that the cohomology type of the two surfaces are ^determined 
by the mapping homotopy type. The zero points are located, the 

30 surfaces are decomposed to patches using gradient lines through 
the zero points. Each of the patches is conformally mapped to a 
rectangle in the parameter domain. To obtain the map between the 
surfaces, these patches on the plane are then matched. 
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Once the mapping between key shapes is known, the points ' on 
the key shapes are selected to serve as the control points . A 
BSpline is used to generate smooth transitions among key shapes. 
This is depicted in Fig. 15 in which a human female face is 
5 morphed using conformal structures into a human male face. In 
this way, we can animate any arbitrary shape. This is especially 
useful for human actors. Facial expressions, gestures, and skin 
deformations of actors at different ages can be stored in a 
database. These stored geometric data can be animated to form 
10 virtual actors. 

Texture mapping without distortion 

Texture mapping of surfaces is very important in both the 
computer gaming industry and the movie industry. The rendering 

15 speed . of a surface is determined by, among other factors, the 
complexity of the geometric model being displayed. For real time 
applications, such as a computer game, simple models are typically 
preferred. In order to improve the perceptual quality of the 
image, images are pasted on the geometric surface using a process. 

20 referred to as texture mapping. 

For a curved surface, texture mapping introduces some 
distortions in the displayed image. The most challenging task for 
introducing texture is to avoid distortion between textures in 
plane and on the curved surface. In industry, geometric modelers 

25 and texture designers are typically different professionals with 
different expertise. Because texture mapping needs to modify both 
geometry and texture, the coordination between these two different 
skill sets are usually difficult and time consuming. 

As discussed above, conformal parametrization has no local 

30 distortion. Using the techniques described above, a geometric 
modeler and a texture designer can integrate their skills easier 
and more efficiently than before. 
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Texture synthesis using Dirichlete Method 

Texture synthesis aims to generate textures to cover a given 
surface from a small texture sample. This is an important 
consideration for graphics design, the movie industry and the 
computer gaming industry. 

Using conformal parameterization, the difficult problem of 
texture synthesis on a geometric surface can be converted into an 
easier problem of texture synthesis on a plane. Using the 
conformal factor analysis and techniques described above, the 
stretching of the texture displayed on the surface can be 
controlled and the geometric properties of the texture on the 
surface can be accurately predicted. 

In order to make the synthesized texture globally smooth, we 
use the Dirichlete method to diffuse the boundaries of texture 
patches. This will make the texture more natural and smoother. 
First the disjoint texture patches on the parameter plane with 
controlled stretching effect are determined. These patches are 
grown until their respective boundaries meet but do not overlap. 
The boundaries of these patches are fixed and the Dirichlet 
problem is solved in the uncovered regions on the surface. Each 
of the color channels are treated as a function, and the 
techniques described above will provide solutions that will 
provide for- a global smooth texture on the surface. 

Volumetric Harmonic Mapping 

Given a 3D manifold, M, a map f z M -> jR 3 is desired that 
minimizes the harmonic energy. In this way the volumetric mapping 
of the original 3D mainfold can be studied in a canonical space. 
The Harmonic energy for f : M-» R 3 is defined as 
E{f)-pffdS M . (70) 

M 

For a discrete system, the harmonic energy is defined as 
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*(/)= (71) 

where k w = ~k^ ff oot(0) , 0 is the dihedral angle opposite to the given 
edge and 1 is the edge length. 

The conjugate gradient method can then be used to minimize 
5 the harmonic energy in order to obtain the harmonic mapping. A 
volumetric harmonic map can be found to map a genus-zero 3D object 
onto a sphere. For the canonical circles on the sphere, a closed 
simple curve on the genus zero object can be found. A Plateau 
problem on the curves can be solved for a conformal deformed 

10 metric. In this way a canonical description of the volume 
enclosed by a surface can be obtained. 

Harmonic mapping is also a useful tool in surgery simulation 
and planning. A physician can construct a 3D brain volumetric 
model from one or more MRI images of the body area of interest. 

15 These MRI images can be mapped onto a 3D sphere. The Physician can 
build a 3D atlas of the body area of interest and compare the 3D 
volumetric data of the new patient's body area of interest with 
the existing atlas data. Because harmonic mapping is unique, this 
technique is a useful method to register brain volumetric data and 

20 would be useful for developing surgery simulations. 

Those of ordinary skill in the art should further appreciate 
that variations to and modifications of the above-described 
methods may be made without departing from the inventive concept 
disclosed herein. Accordingly , the invention should be viewed as 

25 limited solely by the scope and spirit of the appended claims. 
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CLAIMS 

What is claimed is: 

1. A method for matching first and second surfaces, the method 
5 comprising: 

receiving first and second mesh representations of said 
first and second surfaces, respectively; 

conformally mapping both of said first and second mesh 
representations, to a first canonical parameter domain forming 
10 first and second mapped surfaces; 

computing first and second conformal parameter izat ions of 
said first and second mapped surfaces, respectively; 

computing first and second level sets of Gaussian curvature 
and mean curvature for said first and second mapped surfaces, 
15 respectively, wherein said first and second level sets of Gaussian 
curvature and mean curvature are a function of said first and 
second conformal parameterizations, respectively; and 

comparing said first and second level sets of Gaussian 
curvature and mean curvature and in the event that the comparison 
20 exceeds a predetermined threshold, declare a match between said 
first and second surfaces, otherwise declare a mismatch between 
said first and second surfaces. 
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